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Abstract 

It is shown that it is possible for bound fermions on a cosmic string to form a 
superconducting state. Due to the attractive force between them, particles moving 
in opposite directions along the string form bound pairs. This involves a similar 
mechanism to superconductivity in metals at low temperatures. The method of 
Gorkov is used to analyse the system. In contrast to the situation in metals, the 
unusual properties of the string fermion spectrum allow a massless Abelian gauge field 
to provide the required attractive force. This results in far stronger superconductivity 
than usual. A massive gauge field can also be used, in which case the standard results 
apply. 
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1 Introduction 



It is widely believed that during its early evolution the Universe passed through several 
symmetry breaking phase transitions. It is possible that topological defects were formed at 
some of these transitions. The most cosmologically significant of these are cosmic strings [|TJ . 

As well as the strings themselves, it is also possible for particles trapped on the strings 
to produce interesting cosmological effects. A well known example of this is superconduc- 
tivity 0. If an electrically charged scalar field couples to the string Higgs in the right way 
it will gain a non-zero expectation value in the core of the string. Excitations of this con- 
densate act as superconducting currents. The currents are conserved, so the string acts as 
a perfect conductor. In addition to this magnetic flux will be excluded from the string core 
(the Meissner effect). This allows the string to produce irregularities in galactic magnetic 
fields. One criticism of this model is that there are no obvious candidates for the charged 
scalar field. 
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It is also possible to use fermions as the charge carriers. If a fermion field get its mass 
from the string Higgs field, then massless states will exist in the string core 0. As in the 
bosonic case, these fermion currents are conserved, so the string is a perfect conductor. 
However unlike the bosonic currents, there is no Meissner effect, and so the string does not 
act as a superconducting wire in this case (although such currents are sometimes wrongly 
referred to as superconducting in the literature). 

In this paper I will show that it is possible to have a truly superconducting cosmic string 
with Fermi charge carriers. This is achieved by a similar mechanism to superconductivity 
in metals at low temperatures. 

It is a generic feature of an interacting Fermi gas is that if there is an attractive force 
between the particles, the system will have a superconducting ground state ||, |J. In 
this state, particles near the Fermi surface, with opposite momenta and spin will form 
bound pairs, called Cooper pairs 0. Formation of these pairs lowers the energy of the 
system slightly, thus the superconducting state is favoured over normal one. Bardeen, 
Cooper and Schrieffer originally showed the existence of superconducting states with a 
variational argument H. In this paper I will follow the approach of Gorkov M, in which 
superconductivity is analysed using the Green's functions of the theory. 

In a metal the necessary attractive force is provided by phonons. In the cosmic string 
model I will consider it is provided by a massless gauge field. If we were simply dealing 
with a Fermi gas, this type of interaction would not result in a superconductivity. For 
example Quantum Electrodynamics does not have a superconducting ground state. This 
is due to the fact that the photon also produces a repulsive force between like particles. 
However, the unusual properties of the fermion spectrum in a cosmic string background 
alter this fact. 

2 Effective Two Dimensional Action 

Throughout this paper I will denote the time and position on the cosmic string with the 
2- vector x = (t,z), and represent the coordinates in the plane perpendicular to the string 
by x. Similarly k = (k , k z ) gives the energy and momentum along the string. I shall use 
the metric g^ v = diag(l, —1, —1, —1) 

Topologically stable cosmic strings can form at phase transitions whose vacuum man- 
ifolds are not simply connected. The simplest example of this is a U(l) — > I symmetry 
breaking. If $ is the Higgs field responsible for the transition (with ($) = v being the 
usual vacuum solution) and X M is the U(l) gauge field (with coupling g x ), then 



is a cosmic string solution. The masses of $ and X M will be denoted and mx- The two 
profile functions f(r) and a(r) are zero at the string centre (r = 0) and both tend rapidly 
to 1 for r greater than m7 and respectively. 




(1) 
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If any fermions get their masses from $, their spectra will have extra states which 
are bound to the string. I will consider an Abelian string model with two such Dirac 
fermion fields, \l/±, and an additional massless gauge field A^. At finite particle density 
the Lagrangian of the system is 



R 



9- 



— &* A V dp A v + £ s tring 



(2) 



where *f?„' R 



|(1 =F 7 5 )^ s , and £>^f" R = (d„ - ig x Q^ R X^ - igq s A^> R . Gauge invari- 
ance implies Q± — Q± = ±1. The usual masses of the fermions are thus m± = g±v. /i± 
are the chemical potentials of the two fermion species. 

As well as the usual excitations of the fermion fields there are massless states which 
are bound to the string ||. There will also be some massive bound states. The radial 
dependence of the massless states does not depend on their momentum, and so this part 
of the fermion spectrum can be expressed as 



ty s (x,x) = ip s (x)u s (r) + (massive states) 



(3) 



If m± <C 171^ and mx ~ rn<t> then, in the chiral representation of the Dirac matrices, the 
spinors u±(r) can be approximated by 



u+(r) 






V i J 



2iir 



-m+r 



1 
1 

V o J 



2irr 



(4) 



This part of the fermion spectrum has several interesting features. All the massless excita- 
tions of a given field have the same spin and move in the same direction. Which direction 

(with s = ±). 



su f 



depends on how the field couples to <3>. The spinors u s satisfy 7 7 u 
The excitations of the two fields in the above model move in opposite directions and have 
opposite chirality. 

We will be mainly interested in the the massless fermion states, and so for simplicity 
I will ignore the massive part of the spectrum. This will be reasonable if the particle 
numbers and temperature are low enough for there to be no significant occupation of the 
massive states. An effective fermion action can be obtained by integrating out the photon 
field. Due to the form of the states (f|) it is straightforward to remove the dependence 
on the transverse co-ordinates x. The resulting effective two-dimensional fermion action is 
then 

S cS = fd 2 xJ2 4>t [idt + sid z + fi.) ij s +l Yl I d 2 xd 2 x'V s %x-x')ijt(x)iji(x')Mx')Mx) 

J s=± 1 s,a=± J 

(5) 
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where the interaction between the fermions is 
V sa (x — x') = —g 2 q s q a J rdrdOr' dr' d9' u s (r)^u s (r)ua(r')^ v u a (r')D ^(x — x' , x — x ) (6) 

D liV (x — x' , x — x') is the photon propagator. When takes the form (H), the $ dependent 
mass terms in (0) cancel, so they do not feature in (||). 

In a combination of momentum and position space the free-field propagator is 

D S "(k, 3 - g) = r J ^ E ^KM ^^'^.j, (7) 

The dependence on the transverse position has been expressed in polar coordinates (r = |x|) 
using a Bessel function expansion of the plane waves. 

Substituting (|7|) into (§) gives the leading order fermion interaction. This can be 
simplified using 7°7 3 u ± = ±u±. Since Dq U oc g^, the terms in (Bp which correspond 
to interactions between particles of the same species (i.e. when s — a) cancel. Thus 
V sa (k) = (l-8 sa )V (k),with 

V (k) = -g 2 q + q- I k2 _\ 2 _ ie (/ \u + {r)\ 2 J Q (k T r)rdr^ (| \u-{r')\ 2 J Q {k T r')r' dr^j k T dk T 

(8) 

Since the r and r' integrals do not evaluate to delta functions, transverse momentum is 
not conserved in this interaction. This is not surprising since the string breaks Lorentz 
invariance. If q+q- < 0, the force between fermions of different species will be attractive. 
There will be no repulsive force between like particles since they do not interact (at least 
not at tree level). 

An estimate of (H) can be found using the approximate wavefunctions (^). If m± = m, 
it evaluates to 

w) = u {k i - ui) = ^. ^ b (^m^) (8) 

This is the effective potential between the massless fermion states on the string. Non- 
conservation of momentum changes the pole of the propagator to a logarithmic singularity 
at k 2 = k 2 z . For later convenience I will define g = g^J—q^q + /(2n). 

3 Gorkov Functions and Superconductivity 

The attractive force between particles moving in opposite directions with opposite spins in 
the above model suggests that the system will have a superconducting ground state. We 
will be interested in finite temperature as well as zero temperature effects. I will use the 
imaginary-time finite temperature formalism which is easily obtained by defining r — it 
and k = iu n = i(2n+ X)ifT (for fermions), making the replacements (27r) _1 / dk — > TJ2 n 

1 IT 

and J dt — > / dr, and then analytically continuing to real r and integer n. 
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The field equations obtained from ([5]) are then 

\-d T + sid z + Us] ip s (x) = - J d 2 yV(y)4> j _ s (x - y)4>~ s (x - y)ip s {x) (10) 

\-d T + sid z - n.\ tPl(x) = J d 2 yV(y)^ s (x)^_ s (x + y)tp. s (x + y) (11) 

The usual finite temperature Green's functions are defined by 

G s (x — x) = — (T T ip s (x)ipl(x')) (12) 

where T T indicates time ordering with respect to r and (...) indicates the ensemble average. 
It is also necessary to consider anomalous Green's functions [|J, or Gorkov functions 

F s {x - x') = (T T iP_ s {x)iP s {x')) (13) 

F}(x-x')^-(T T ^l s (x)4(x')) (14) 

In the vacuum, or a state with a definite number of particles, F will be zero. More generally 
this need not be the case. 

Using ([Top and (|TT|), and the finite temperature version of Wick's theorem |5| 

[-d T + sid z + Vs]G s (x) + J d 2 yV(y) [G- s (Q-)G a (x) + F s (y + )F}(x + y)] = 5(x) (15) 

[_ 9r _ Sldz _ n_ s ]Fl{x) - J d 2 yV(y) {G- s (0-)F}(x) - F}(y + )G s (x - y)] = (16) 

here 0^ denotes the 2-vector lim r? ^ (±?7, 0) with rj > 0, and y ± denotes y + ± . If the 
average number of ip s particles per unit length is N s , then (ip\ip s ) = G s (0~) = N s . 

Changing to Fourier space, with G(k) = J dzdre luJnT ~' lkzZ G(x), and similar expressions 
for F\ and V, 

[iu n -e k )G s (k)+A s (k)F}(k) = l (17) 

[iu; n + C s k ]FKk)+^s(k)G s (k) = (18) 

where ^| = sk z + S s — /i s . The self energy of the fields is S s = —VN_ S , where 
V = V(k = 0) =Jd 2 xV(x). 

Expressing G(0~) in momentum space gives 

r Tdh 

N s = J i^Y.^ nV Gs{k) (19) 

which can be used to determine the chemical potentials, fi s . The rj — > + limit is implied in 
this and the following expressions. As it stands, the above expression for N is not correct, 
since the vacuum contributions need to be subtracted off. This is obtained by evaluating 
( fOP when jj, and T are zero. 
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The gap function A s (k) is defined by 

A.(*) = / ^ E ^^(p - WCp) (20) 

It will be zero in a non-superconducting state. 

For simplicity I will consider a system with equal numbers of left and right moving 
particles (N s = N). In this case fi s = /i, the Green's functions of the two particle species 
are related by G s (k) = G(k , sk z ) and F s {k) = sF(ko, sk z ). Similarly A s (k) = sA(ko, sk z ). 
The equations (|T7D and are then solved by 

Gik) = -^WF (21) 



*W = , fxTTXT^ (22) 



A(*Q 
w* + # + |A(A0| 2 

If A 7^ then it can be seen from that the there is an energy gap in the fermion 
spectrum at the Fermi surface = 0). This is not the same as a mass, since it is not at 
the zero energy part of the spectrum. The Green's function (|2l| ) gives the excitations of a 
state with an average of N particles of each species. If A ^ these are not particles and 
holes (as in the normal ground state), but combinations of them, called quasiparticles. 

Substituting (|22| ) into (|20| ) gives an expression for A(k). The only physical quantities 
are those involving on-shell values of k, which can be found by analytic continuation to 
real values of iu n = ko- For the filled states the on-shell values of fco are — Ek with E^ = 



£l + \A k \ 2 and A k = A(-E k , k z ). The empty states have k = E' k with E' k = ^ + |A' fc | 2 
and A' fc = A{E' k ,k z ). 

The gap equation (f20D for A^ is then 



A k = f ^IL^n(-E p )U[2E p E k - - |A P | 2 - |A fc | 2 ] 

~ ^Fri(E' p )U[-2E' p E k - 2& k - \A'/ - |A fe | 2 ] (23) 



where n(E) = [exp(E/T) + l] -1 . U is defined by (||). The corresponding equation for 
A^ can be obtained with the replacements A^ — > A^. and E k — > —E' k . In the absence of 
magnetic fields, we can take A& to be real and positive. 

Clearly A(k) = is a solution of (|23|). This corresponds to the normal ground state of a 
Fermi system. At zero temperature this solution has the fermion states (completely) filled, 
in order of increasing energy up to the Fermi surface. It is also possible that a A(k) ^ 0, 
superconducting solution also exists. In this case (at T = 0) the fermion states either side 
of the Fermi surface are only partially filled. Thus the Fermi surface is no longer sharp in 
the superconducting state [||]. 
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4 The Gap Equation at Zero Temperature 



At zero temperature fl23|) simplifies, and it is possible to remove the k dependence from 
the integral. The gap equation can then be solved by constant A and A'. At T = 0, 
n(E) = 9{—E) where 9 is the step function, so the second term of (|23| ) vanishes. The first 
can be evaluated using the change of variables = A sinh 2k and £ p = A sinh 2\- The gap 
equation (|23| ) then reduces to 

A = A / ^U[AA 2 sinh 2 (x - «)] (24) 

J 27T 

Similarly, using = — |A'| sinh 2k, the corresponding equation for A' becomes 

A' = A [ ^C/[-(|A'| - A) 2 - 4 1 A' | A cosh 2 (x - «)] (25) 

J 27T 

The dependence on k (and hence /c) is then simple to remove. Unfortunately analytic 
evaluation of the resulting integrals does not seem possible. However we expect A to be 
very small. This fact can be used to approximate the integrands of (0,^) • For simplicity I 
will take m + = m_ = m, in which case U[—(p — k) 2 ] is given by (0). The main contribution 
to the integrals occurs for moderate values of x, when e~' x ' <C e'^ <C 2m/A. In this region 
the denominator of @ is approximately constant and 2 sinh x ~ e' x L For \x\ > ln(2m/A) 
the integrand of (|24]) tends to zero rapidly. This suggests the approximation 



;If/[(2Asinh X ) 2 ]^ 2 

Z7T 



2m 

l n ~n X 

A 



e(ln^-|x|) (26) 



Using this approximation (24) becomes 



/ .in(2m/A) / 2m , A , A , / 2m\ 2 . . 

A^Ag 2 i n x dx = A/ In— 27 

J-ln(2m/A) V A / V A / 



Thus either A = (normal solution) or 



A re 2mexp ( — | = 2mexp I -J ^ -) (28) 

which is the superconducting solution. This value of A contrasts with the more usual result 
of A ~ exp(— 1/g 2 )- It can be seen that ( p8|) has an essential singularity at g = 0. This 
is a generic feature of superconductors, and explains why they cannot be analysed with 
standard perturbation theory. 

To leading order the approximation (|26|) also works for (^), and so A' w A. There 
will be a small phase difference between the two constants because the right-hand side of 
(^5|) includes a small imaginary part. 

The total energy of the state is given by 

E= lim X>'0* + £jG,(ar) (29) 
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When there are equal numbers of ip+ and ip- particles, the momentum space version of 
( p9|) reduces to 

Tdk 2 



(30) 



where G(k) is given by (pl|). Using (|19"1), the chemical potential is found to be [l = 
2ttN — VN. As it stands, for A^O, (^) is divergent. If we introduce a cutoff at = A 
in the divergent part of the integral, it evaluates to 



E = 2ttN 2 



2VN - — A 2 

47T 



1 + ln 



4A 2 
A^ 



(31) 



In the calculations of section [5] all energies and momenta were assumed to contribute 
to the superconducting condensate. In practice contributions from the massive part of the 
fermion spectrum (which were ignored for simplicity) will alter the effective potential for 
> m, producing some sort of cutoff, probably around the off-string fermion mass. This 
suggests that m is a reasonable value for A. In this case 



E = 2irN 2 - 2VN A 2 

47T 



2 

1 + - 

9 



(32) 



Clearly (whatever the value of A), fl32|) is smaller when A 0, thus the superconducting 
ground state is favoured over the normal A = one. 

Superconducting cosmic string models have been studied previously using a fundamen- 
tal bosonic scalar field as a superconducting condensate M. A connection with this work 
can be made by considering an effective bosonic field ip+ip- (not to be confused with ipip^). 
This will have a non-zero expectation value in the superconducting state. This is equal 
to F(x — > 0). Evaluating with ( p^) gives a divergent integral. As with ([32]), it is 

reasonable to impose a cutoff. In this case (ip-4>+) ~ A/g, which is very small. 



5 Critical Temperature 

At high temperatures, thermal excitations prevent the fermions forming the bound pairs 
required for the superconducting state. As the temperature is increased from zero, a 
progressively smaller fraction of the fermions will be paired up. At some critical tem- 
perature T c there will be no Cooper pairs left and the system will just be a normal (i.e 
non-superconducting) gas at finite temperature. 

Solving ( p3| ) when T ^ is much harder than the zero temperature case. It is no longer 
possible to remove the k dependence from the integral, and A is no longer a constant. 
Despite this it is still possible to estimate the critical temperature, T c , above which the 
ground state of the system is no longer superconducting. This can be found by finding the 
value of T which satisfies ( |2"3"D when A — > 0. 

While taking A small simplifies (|23|), the k dependence of the integral still creates 
problems. To estimate T c I will take A& and A^ to be equal and constant. Using the 
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approximations in the previous section, the gap equation when = reduces to 

p\n(2m/A) / 2m \ [ A i\ 

\^f\ d X In— Uanh — e W (33) 

J-ln(2m/A) V A / \ll c J 

where £ p = Asinh2x. Changing variables to w = A/(2T c )e x , (|33|) becomes 



m/Tc , f T c \ tanhw , , , 

, In — w dw 34 

g z Jo \m ) w 

The leading order (in T c /m) contribution of the integral can be found by integration by 
parts 

1 /, m" 



9 2 V T c 



-In- (35) 



Solving this gives 



T c ~ mexp f — ) ~ A| T=0 (36) 
V 9) 

For £k 7^ there will be additional terms of order In m/T c on the right hand side of fl3"5]). 



Since m 3> T c , these can be neglected to first order and so ignoring the k dependence of 
A is a reasonable approximation. 

The corresponding result for metals is also T c ~ A, although in this case A is far smaller. 
Since T c m, the cosmic string will not become superconducting until the universe has 
cooled considerably. The temperature at which this occurs is extremely sensitive to g. For 
example, suppose that \I/± are GUT scale fermions with m ~ 10 16 GeV. If g = 1/10 then 
T c ~ 10 n GeV, while if ~g = 1/100, T c ~ 10~ 28 GeV. It is therefore possible to construct 
a model which becomes superconducting at a temperature far below the string formation 
scale, without having to introduce unnaturally small fundamental constants. Note that it is 
still possible to get the standard non-superconducting fermion currents at any temperature. 
All that is required is different numbers of left and right moving particles on the string. 



6 Higher Order Effects 

So far I have ignored higher order corrections to the interaction. These will be particu- 
larly significant for calculating the self energy of the fermions. To leading order this is 
proportional to V = V + ~{k = 0), which according to (|9|) is divergent. 

Taking into account loop corrections from the fermion states on the string the potential 
V is given by 

V SC = V SC + ]T V sa U ab V bc (37) 

a,b=± 

where V Q SC is the free field potential ([|) and U sr is the polarisation. To leading order in g 

72 

IU(£;) = -5 ab [ j^G a (k + p)G a (p) (38) 
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At zero temperature the contribution from the massive fermion states will be the same as 
in the vacuum. Since (after renormalisation) these do not contribute to V when k = 
they can be ignored. 



To leading order in g, we can take A = in (|38). In this case the Green's functions 
reduce to G(k) = {iuj n — ^k) 1 - There are also contributions to fl38|) involving F and F\ but 
since they are 0(A 2 ) they can also be dropped. With these approximations (j3~8| ) evaluates 
to 

U<*(k) = ^5 ab (39) 
at any temperature. Substituting this into ( p7|) and solving gives 

W -vr - v„ 2 {W> 

Thus V + ~ — > as k — > 0, so the terms considered in section |3] do not contribute to the self 
energy of the fermion fields at all. 

The corrections to V also create a pole at k 2 z — k\ = 4m 2 e -1 ^ . The residue is of order 
g _2 e _1// ^ <C A, and so the pole does not significantly alter the results of the previous 
sections. 

In addition it can be seen that the interaction between fermions of the same species is 
no longer zero. This will lead to further corrections to the self energy of the fermions. To 
fully analyse this, Green's functions of the form 

G s> - S (x - x') = -{T T ^ s {x)^_ s {x')) (41) 



need to be considered in addition to the usual Green functions ([12]). The main contribution 
to G s - S will come from particle-antiparticle interactions with energies near to zero, which 
may give the Green's functions an effective mass term. The superconducting effects result 
from quasiparticle interactions with energies near the Fermi surface i.e. in a different part 
of the spectrum. The corrections arising from (|4l|) are therefore unlikely to affect the 
superconductivity of the system and I will not consider them further in this paper. 

It is also possible that the photon could gain a mass, perhaps from some other fields 
not bound to the string, or from the massive \& states when T ^ 0. In this case a mass 
term will appear in (^). Alternatively the the force could be provided by a massive boson, 
possibly the string gauge field X. If A were used there would be contributions to the 
potential from string bound states as well as the continuum states. For a constant mass, 
the potential will resemble (H) or @, with k\ replaced by k\ + m\. In this case the self 
energy of the particles is 

E = -^f-^-lV 1 ln^ (42) 
Y Am 2 J 2m 

If rriA 3> A, fl26| ) is not a good approximation of V and 

±U[m\ + (2Asinh X ) 2 ] - f ^ (in ^ - | X |) (43) 

V 2m J 1 
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is more appropriate. In this case the non-zero solution of the gap equation is 



A 



rriA exp 



'"A 
2 m 



"I 1 



V 



In (it) 



2 ~9 2 



(44) 



which resembles the more usual result for metals. Strictly speaking, if A is massive the 
cosmic string is no longer superconducting, since the fermions no longer carry (unbroken) 
charge. However \l/± could carry other charges in which case the system will superconduct. 
Of course this assumes that interaction with the other gauge bosons is not strong enough 
to cancel the attractive force arising from the field A. 



7 Conclusions 

It has been shown that it is possible for a cosmic string with massless fermion bound states 
to have a superconducting ground state. In order for this to happen there must be fermions 
moving in both directions along the string, with an attractive force between them. 

If the oppositely moving fermions have opposite charges, the attractive force can be 
provided by a gauge boson with zero mass. This is in contrast to an interacting Fermi gas 
in a spatially homogeneous (i.e. no cosmic string) background. Metals at low temperature 
fall into this category. In this case a massless intermediate boson will not give a super- 
conducting ground state. The unusual features of the string fermion spectrum remove the 
repulsive part of the interaction, which would normally prevent superconductivity. 

One characteristic of a superconducting state is an energy gap A in the excitations of 
the state. The occurs at the Fermi surface, and is usually of order exp(— l/g 2 ), where g is 
coupling strength of the gauge boson providing the force between the fermions. If the gauge 
boson is massless and the fermions are bound to a cosmic string, A is far higher, of order 
exp(—l/g). This is lower than might be expected, since non-conservation of momentum 
by the string reduces the force. 

The quantity A also determines the critical temperature of the system, above which 
it ceases to be superconducting. Since A is very small, this will be far lower than the 
fundamental energy scale of the model. 

The electromagnetic effects of the superconductor are proportional to the number of 
superconducting fermions in the system. At zero temperature this is all of them, so the 
electromagnetic effects need not be small even though A is. These can be analysed by, for 
example, perturbing the Green's functions [f|] (if the field is weak). 

Some possible candidates for the fermions to make up the condensate are superheavy 
GUT fields, such as those occurring in E 6 or SO(10). 

Eq has superheavy versions of electrons, neutrinos and down quarks. An attractive 
force could be provided by GUT interactions. Since these would involve massive bosons, 
A would be given by ( fPf ) rather than (|25|). Alternatively if the force were provided by 
ly-bosons, the massless boson result fl28|) would apply, at least until the electroweak phase 
transition. 
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The only fermions coupling to the string Higgs field in an SO(10) GUT are right-handed 
neutrinos. Since at least two species are required to form Cooper pairs, this type of GUT 
string does not have superconducting states. 

Interestingly, electroweak strings have u, d and e massless string bound states. Conve- 
niently the up quarks (positively charged) move in the opposite direction to the electrons 
and down quarks (negatively charged), so these strings have all the required features of 
the model discussed in this paper. Unfortunately electroweak strings do not appear to be 
stable. However the effects of the fermion fields have not been fully investigated in these 
models, so they could be stable under the right conditions. 

Acknowledgements 

I wish to thank Warren Perkins and Simon Hands for useful discussions, and PPARC for 
financial support. 

References 

[1] A. Vilenkin and E. P. S. Shellard, Cosmic Strings and other Topological Defects (Cam- 
bridge University Press, Cambridge, England, 1994) 

[2] E. Witten, Nucl. Phys. B 249, 557, (1985). 

[3] R. Jackiw and P. Rossi, Nucl. Phys. B 190[FS3], 681, (1981). 

[4] A. A. Abrikosov, L. P. Gorkov and I. E. Dzyaloshinski Methods of Quantum Field 
Theory in Statistical Physics (Dover , New York, 1963). 

[5] A. L. Fetter and J. D. Walecka Quantum Theory of Many-Particle Systems (McGraw- 
Hill, New York, 1971). 

[6] L. N. Cooper, Phys. Rev. 104, 1189, (1956). 

[7] J. Bardeen, L. N. Cooper and J. R. Schrieffer, Phys. Rev. 108, 1175, (1957). 

[8] L. P. Gorkov, J. Exptl. Theoret. Phys. 34, 505 (1958) [Sov. Phys. JETP 34, 735 
(1958)]. 



12 



